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Abstract: This paper presents a compression method for inertial and magnetic sensor data, 
where the compressed data are used to estimate some states. When sensor data are bounded, 
the proposed compression method guarantees that the compression error is smaller than a 
prescribed bound. The manner in which this error bound affects the bit rate and the estimation 
error is investigated. Through the simulation, it is shown that the estimation error is improved 
by 18.81% over a test set of 12 cases compared with a filter that does not use the compression 
error bound. 
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1. Introduction 

Largely because of the MEMS technology, inertial sensors (accelerometers and gyroscopes) are be- 
coming smaller and cheaper [1], which makes it possible to use inertial sensors in many applications. 
Inertial sensors are used in motion trackers [2], personal navigation systems [3] and remote control 
systems [4]. 

In some applications such as body motion trackers (for example, the product 'Moven' by XSENS), 
inertial sensors are used to track body movement. As the number of inertial sensors increases, the size 
of sensor data increases accordingly. The sensor data is transmitted to the microprocessor board through 
wired or wireless communication channels. In a wireless communication channel, the transmission speed 
is relatively low compared with a wired communication channel. The size of the sensor data needs to be 
reduced if it exceeds the capacity of the communication channel. For example, the transmission rate of 
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raw sensor data for one MTx (commercial inertial and magnetic sensor) could be as high as 72 Kbps: 9 
sensor outputs (3 accelerometers, 3 gyroscopes, and 3 magnetic sensors) x 500 Hz (maximum sampling 
rate) x 16 bits (16 bit A/D conversion for each sensor). If four MTx are used, the size of the sensor 
data exceeds the capacity of Zigbee (maximum rate is 250 Kbps [5]). More applications are expected as 
networked control and monitoring systems are becoming more popular [6]. 

One way to reduce the size of the sensor data is to compress the sensor data before transmission and 
decompress the received data in the microprocessor board. Data compression has been extensively stud- 
ied in many areas [7]. In applications such as body motion trackers, real-time compression is preferable, 
in order to avoid delayed sensor data transmission and consequently the delay in motion estimation. 
One of the most popular real-time compression methods is ADPCM (Adaptive Differential Pulse Coded 
Modulation) [8], which is optimized for voice data. In [9], a simplified ADPCM method is used for 
inertial sensor data compression, where the maximum error (the difference between the original data and 
the compressed-and-then-decompressed data) is only relatively bounded (e.g., 1% of the sensor data). 

The performance indices of data compression are the bit rate and the quality of compressed data. In 
voice data compression, the quality of compressed data is evaluated by listening to the compressed-and- 
then-decompressed voice data. This rather subjective evaluation makes sense since the final destination 
of compressed data is the human ear. On the other hand, the final destination of compressed inertial 
sensor data is usually a filter (such as a Kalman filter), where orientation is estimated. Thus the quality 
of compression should be judged by how the compression affects the estimation error. 

In this paper, a modified ADPCM method is proposed, where the absolute maximum error bound is 
explicitly given. Also, we investigate how this error affects the estimation error. A part of this paper was 
presented in [10]. 

2. Inertial and Magnetic Sensor Data Compression and Estimation 

The overall process of compression and estimation is given in Figure 1, where A; is a discrete time 
index. The objective is to estimate some states x{k) (attitude, heading, position, etc.) using inertial 
sensor data y{k) at a limited data transmission rate. The inertial sensor data y{k) is compressed into 
d{k) and transmitted to the microprocessor board. The compressed data is decompressed into y{k) and 
the state x{k) is estimated using a filter. 



Figure 1. Overview of inertial and magnetic sensor data compression and estimation, 
data transmission 
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Since the objective is to find a good estimator of x{k), the quality of compression is considered good 
if the estimation error x{k) — x{k) is small. The quality of the compression algorithm is evaluated using 
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the following estimation error covariance: 

Fervor = ^{{X - X{y))' {X - X{y))] (1) 

where x{y) is an estimator when y is used as an output. 

Note that Perror depends on the filter algorithm used to compute xiy) in addition to the compression 
algorithm. 

The ideal compression algorithm minimizes both Perror and the bit rate. However, usually if Perror 
is small, the bit rate tends to be large. In Section 4, we propose a compression algorithm where the 
maximum compression error is bounded. The maximum compression error bound plays a role of design 
parameter to adjust Perror and the bit rate. 

3. Modified ADPCM Algorithm 

The ADPCM block scheme is given in Figure 2. We assume that y{k) is the output of Uy bit uniform 
quantizer, where y{k;) satisfies 

(2) 



Figure 2. Encoder and decoder block scheme, 
compression decompression 




Let the quantization size 5 of y{k) be defined by 



6 



Vn 



2"!/- 



(3) 



If there is more than one sensor, we need one encoder for each sensor. 

The sensor data y{k) is compared with the predictor output y{k). The difference d(k) is coded into 
d{k) and this d{k) is transmitted to the estimator board. In the standard ADPCM, d{k) is a quantization 
index i{k). In this paper, d{k) consists of one bit mode information m(k) and a quantization index i(k): 



d{k) 



m{k) 

m 



(4) 
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In the decoder, the decompressed data is y{k) = y{k)+d{k). The predictor output y{k) can be computed 
from d{k) and thus does not need to be transmitted. 

The adaptive predictor uses the same pole-zero configuration as that in CCITT G.726 ADPCM, which 
is an ADPCM speech compressor/decompressor protocol proposed in 1990 [11]: 

2 6 

y{k) = Y,a,{k- i)y{k -i) + J2b,{k- l)d{k - z) (5) 

1=1 i=l 

From the assumption (2), if y{k) > y^ax from (5), we set y{k) = y^ax- Similarly, if y{k) < -y^ax, we 
sot y{k) = -ymax- 

The adaptive algorithm in the G.726 protocol is used to adjust Oj and bi and the detail is given in [1 1]; 
the tone and transition detector part was omitted since the part is only for voice data. 

The compression error ec{k) is the difference between the original signal y{k) and the decompressed 
signal y{k): 

ec{k) = y{k)-y{k) 

= {y{k)+d{k))-{m+dm (6) 
= d{k)-d{k) 

Standard ADPCM algorithms [8] will be modified so that the maximum error is bounded as follows: 

(7) 

Now d(k) coding is explained. The mode bit m{k) in d(k) is used to ensure (7). If the compression 
error edk) of a standard ADPCM method satisfies (7), then the mode is 0 (i.e., m{k) = 0). As will be 
seen in Section 3.1, this is true if \d{k) \ < 2^=('')5, where ysik) is an adaptive scaling factor. On the other 
hand, if the compression error edk) of a standard ADPCM method does not satisfy (7), then the mode 
is 1 (i.e., m{k) = 1) and a special uniform quantizer is used as in Section 3.2. Thus the mode m{k) is 
given by 

r 0 if < 1 

I 1, otherwise 

where ys{k) is an adaptive scaling factor. 

The quantization index i{k) is defined differently when m(k) = 0 and when m{k) = 1. 

3.1. Quantization index when m{k) = 0 

If m{k) = 0, a signal d{k) is quantized with Ud bits with a logarithm quantizer with an adaptive 
scaling factor ysik), where the quantized index i{k) (1 < \i\ < 2"'*^^) satisfies 

The sign of the index i{k) is the same as that of d{k). If d{k) = 0, then i = 1. 

Coefficients fi in (9) are computed from /i law [8] so that fi satisfies the following: 

log,(l + ^ i 

loge(l + /i) 2"<*-l 
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Table 1. fi values for rid = 5. 



/o 


/i 


/2 


/3 


h 


0 


0.0237 


0.0503 


0.0799 


0.1131 
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h 


/r 


h 


h 


0.1501 


0.1915 


0.2379 


0.2899 


0.3479 


/lO 


hi 




/l3 


/l4 


0.4129 


0.4855 


0.5667 


0.6575 


0.7593 


/l5 


/l6 








0.8729 


1.0000 









In this paper, = 5 is used and fi values for the case of = 5 is given in Table 1: 
The scaling adaptation factor ys{k) is computed similarly with the standard ADPCM algorithm except 
that ys{k) is bounded as follows: 

3 < ys{k) < Vs (10) 

We note that yg is chosen so that (7) is satisfied. First we are going to derive the upper bound of ec{k) 
when ys is given. 

The decompressed signal d{k) is computed as follows: 

-fi-l + fi 



where 



The error edk) is then 



d{k) = sgn{t{k))2y'^''^5- 



sgn(a) 



\ec{k)\ = \d{k)-d{k)\ 




(11) 



2 

From the fact that fi is monotonically increasing and (10), we have 



< 2^^ (5 



2 

Given Cmax, to satisfy (7), ys should satisfy the following 

Thus if ys is chosen to satisfy (12), the quantization error is always smaller than Cmax when m(k) = 0. 
We also note that in addition to the global bound Cmax, if index i{k) is known, we have a less conservative 
bound given in (1 1): 

edk) = 2^-W^ ^' (13) 
This bound will be used later in the estimation problem. 
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3.2. Quantization index when m{k) 



From (8), m{k) = 1 if \d{k)\ > 2^"*^'^^. If m{k) = 1, then the logarithm quantizer used cannot 
guarantee the maximum error (7). Noting that d{k) = y{k) — y{k), we can see that the mode is 1 if the 
difference between the output y(k) and the predicted value y{k) is large, which happens when the signal 
change is not smooth but instead rather abrupt. 

To ensure the maximum error condition (7), we introduce a uniform quantizer when the mode is 1. 
An example is given in Figure 3, where y{k) is outside the [y{k) — 52y'''^^\y{k) + 52^"^'^)] interval and 
the mode is 1 . The upper and lower intervals of the mode 1 interval (the mode 0 interval is the shaded 
area) are quantized with a uniform quantizer (quantization size is 2emax)- The formal definition of the 
index in mode 1 is given as follows. Let uiength and kength be defined by 

Ulength = Vmax " {y{.k) + 522^^('=)) 
kength = {y{k) - 52^'^^'^) + y^ax 

Let Mieuei (the number of quantization levels for the upper interval uiength) be defined by 



0, 



if Ulength < 0 



ceil (^^1^^) , otherwise 
ceil(Q;) is the smallest integer no smaller than a. 



where 



Figure 3. If \y{k) - y{k)\ > 2?^=(^)(5, then m{k) = 1. 



y(k) is outside the 
shaded interval 

y{k) ^ 



— y(k) + S2''"' -L 

— m) 

— y{k) - ^2'-<« 

^max 



uniform quantizer 
is used for two 
intervals 



The index i[k) is given by 



floor M iiy{k) > + 52^^W 

^ ^ ^ Ulevel + noorl y-'%:Z7''^ ) ify(fc)<y(A:)-^2^-W ^ ^ 

where floor(a) is the largest integer no larger than a. 
The decomposed signal d{k) is computed as follows. 



d{k) 



evel 



max{ymax, y{k) + 52f=('') + emax + i{k)2eraax}, if i < Ul 

mm{-y^ax, y{k) - 52J'=(^) - e^ax - {i{k) - Uievei)'2emax} , otherwise. 



Since a uniform quantizer is used, the compression error edk) in mode 1 is bounded by 

|ec(/c)| < e^k) = emax (15) 
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Note that the number of bits for the index i{k) is given by 

ceil (log, (ceil (^)+ ceil fl^))) 

When m{k) = 1, ni{k) changes depending on y{k) and y{k). Note that when m(k) = 0, n^^ (the number 
of bits for i{k)) is constant. Also note that rij can be computed in the estimation board and thus does not 
need to be transmitted. 

4. Kalman Filter Compensating the Compression Error 

In Section 3, a compression method is proposed, where the maximum compression error is Cmax- Also 
if m{k) and i{k) are known, bounds of the compression error are given by (13) and (15). In this section, 
we use this information in a Kalman filter. 

We assume that y{k) is generated by a linear system 

x{k + l) = Ax{k)+w{k) 
y{k) = Cx{k)+v{k) 

where x G -R" is the state, y E is the output, and w{k) and v{k) are uncorrelated, zero-mean, white 
Gaussian noises that satisfy 

E{w{k)w'{k)} = Q, E{v{k)v'{k)} = R 

In the standard Kalman filter, x{k) is estimated using y{k). If y{k) is compressed, y{k) = y{k) — ec{k) 
is used instead. From (13) and (15), ec{k) can be computed, which is used to reduce the estimation error. 
If we assume ec,i{k) (z-th element of ec(/c)) has a uniform distribution, E{el,-{k)} = lel^{k). By 
treating the compression error edk) as measurement noise in y{k), the following model can be used for 
an estimator. 

x{k + l) = Ax{k)+w{k) 
y{k) = Cx{k)+v{k) 



(16) 



(17) 



where 



R{k) = E{v{k)v{k)} = R + 



(18) 



Since the compression error is compensated in the estimation algorithm, we can expect a smaller 
estimation error, which is verified through the simulations in Section 5. When a Kalman filter is used for 
(17), the estimation error covariance P{k) = E{(x{k) — x{k)){x{k) — x^k))'} can be computed from a 
Riccati equation [12]. 



P{k + 1) = AP{k)A' + Q- AP{k)C\CP{k)C' + R{k))-^CP{k)A' 



(19) 



Since R{k) depends on y{k), we cannot compute P{k) before simulation. To evaluate the estimation 
error covariance without simulation, we use an upper bound of R{k): 

o2 



R{k) <R + 



-^max,p 



Rr. 
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Using this Rmax in place of R{k), we can find a steady-state solution to (19). 

P = APA' + Q- APC'iCPC + RmaxY^CPA' (20) 

P can be considered as an upper bound of P{k) in (19). From P, we can see how Cmax^i affects the 
estimation error. 

A similar idea is used in [13], where a networked estimation problem is considered. 
5. Simulation 

We compared three data sets using the proposed compression algorithm. Original data is 1600 bits/s 
for each sensor : 16 bit A/D converted data (i.e., Uy = 16) with the sampling rate being 100 Hz. We used 

= 5: that is, the number of bits for the quantization index when m(k) = 0 is 5. Cmax for each sensor 
is chosen so that e^ax = 3005, where 5 is different for each type of sensors, yg = 12 is found to satisfy 
(12). 

Bit rates for the three compressed data sets are given in Table 2. All three data sets are obtained using 
XSENS MTx (3 accelerometer, 3 gyroscopes, and 3 magnetic sensors). Holding MTx with a hand, we 
moved MTx slowly (data set 1) and fast (data set 2). Data set 3 is obtained from a personal navigation 
system, where MTx is attached on the shoe of a pedestrian [3]. 

The bit rates of data set 3 is the largest because the change of data is the most abrupt. In particular, 
when the shoe contacts the floor, there is a large change in the accelerometer and gyroscope outputs, and 
consequently the compression algorithm becomes m{k) = 1 more often. 



Table 2. Bit rates for 3 inertial and magnetic sensor data sets. 





accelerometers 


gyroscopes 


magnetic sensors 




643.2 


641.6 


617.6 


data set 1 


648.0 


632.0 


622.4 




656.0 


628.8 


614.4 




659.2 


660.8 


617.6 


data set 2 


662.4 


640.0 


622.4 




668.8 


638.4 


619.2 




771,2 


715.2 


638.4 


data set 3 


798.4 


696.0 


633.6 




694.4 


686.4 


627.2 



To see how the compression error affects the estimation error, a simple one dimensional attitude 
estimation problem is considered. An attitude (9) is estimated using two outputs: yi is an inclinometer 
output and yg is a gyroscope output, where 



yi = 0 + Vi 

Vg = ^ + ^9 



(21) 
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where Vi and Vg are measurement noises and qi = E{vf} = 0.13 x 10^^ and qg = E{Vg} = 0.78 x 10^^. 

An indirect Kalman filter ([12]) is used to estimate 9 using yi and yg. In the indirect Kalman filter, the 
gyroscope output is integrated to compute 9i, i.e., 

Ot = Vg 

Defining 9s as 

9s = 9i — 9 

we obtain 

9s = Vg (22) 

In the indirect Kalman filter, 9s is first estimated and 9 is obtained indirectly from 9 = 9i — 9s- 
By discretizing (22) with the sampling period T (T = O.Olsec), we have 

9s{k + l) = 9sik) + VTvg ^23) 

di-Va = ds{k) - Vi 

Now the proposed compression algorithm is used to compress y^ and yg. The maximum compression 
errors emax,i (emax for yi) and emax,g (emax for yg) affect both the bit rate and the estimation error. If Cmax 
is small, the chance that the mode becomes 1 increases. Thus the bit rate becomes large. How the bit 
rate changes with the changing Cmax is given in Figure 4. The data yi and yg are generated using Matlab. 



Figure 4. Relationship between Cmax and the bit rate. 




e for inclinometer output 



From Figure 4, we can see that the bit rate of the inclinometer outputs increases rapidly as emax,i is 
decreased. On the other hand, the bit rate of the gyroscope outputs does not change much as emax,g is 
decreased. The bit rate depends on how often the mode becomes 1 : note that rii (the number of bits 
needed for the quantized data when the mode is 1) is generally larger than 12^ (the number of bits needed 
when the mode is 0). If the original signal is sufficiently smooth, d{k) is small since the predicted value 
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y{k) is very close to y{k). Thus even if we decrease Cmax, d(k) still satisfies \d{k)\ < 2^=(''')5 condition 
in (8). 

The yi and yg signals and the compressed signals are given in Figures 5 and 6, where emax,i = 0.2876 
and emax,g = 0.0122. We can see that the gyroscope output is relatively smooth compared with the 
inclinometer output. Thus the bit rate of the gyroscope output is relatively insensitive to the changes 

m ^max,g- 

Figure 5. Inclinometer output, decompressed data, and error. 



„ 0.4 




time (sec) 



Figure 6. Gyroscope output, decompressed data, and error. 
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The effects of changes in e-max on the estimation error are given in Figure 7, where the estimation 
error is predicted using (20). Actual estimation error from simulation is given in Figure 8, where 
\/E{dik) -e{k)y is computed. 
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Figure 7. Predicted estimation error (P) of the proposed filter. 




Figure 8. Estimation error (y I](6'(A;) — 6{k)Y ) of the proposed filter. 




The relationship between bit rates and estimation error is presented in Figure 9, where data are from 
Figures 4, 7 and 8. In the left graph, the points 1-6 have similar bit rates but different Perror- Thus in the 
following simulation, we chose the point 1, which corresponds to emax,i = 0.2876 and emax,g = 0.0122. 
We compared three different filters: (a) a standard Kalman filter using yi and y^; (b) a Kalman filter 
using yi and yg with compression error compensation (proposed in Section 4 ); (c) a Kalman filter using 
{ji and ijg without compression error compensation. We randomly generated 12 data sets and the results 
are given in Table 3. 
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Figure 9. Bit rates and Perror'- predicted value by (21) and the experiment result. 
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Table 3. Bit rates and estimation error of 3 filters: (a) a standard Kalman filter with uncom- 
pressed data, (b) the proposed method with yi and yg, and (c) a Kalman filter with yi 
and yg. 



experiment 


Bit rate 


p 

^ error 


% improvement 
((c) -(b)) /(c) 


Va 


m 


(a) 


(b) 


(c) 


1 


684.8 


620.3 


0.318 


0.437 


0.452 


3.32 


2 


654.4 


623.3 


0.173 


0.361 


0.492 


26.62 


3 


640.0 


621.3 


0.138 


0.345 


0.402 


14.19 


4 


645.7 


622.2 


0.139 


0.312 


0.455 


31.39 


5 


631.2 


625.7 


0.161 


0.382 


0.447 


14.50 


6 


633.6 


621.6 


0.121 


0.310 


0.422 


26.55 


7 


692.7 


692.7 


0.308 


0.427 


0.445 


3.96 


8 


640.8 


622.2 


0.120 


0.280 


0.317 


11.77 


9 


634.4 


618.6 


0.134 


0.295 


0.537 


45.15 


10 


696.1 


618.8 


0.423 


0.487 


0.503 


3.14 


11 


638.4 


632.9 


0.191 


0.394 


0.464 


14.99 


12 


638.4 


621.5 


0.136 


0.292 


0.417 


30.13 


average 


652.5 


628.4 


0.197 


0.360 


0.446 


18.81 



It is not surprising that the Perror of the standard Kalman filter is the smallest because the original data 
Hi and Ug are used for measurements. We can also see that the Perror of the proposed filter is smaller than 
that of the filter (c). On average, the estimation error of the proposed filter is smaller by 18.81% compared 
with that of the filter (c). Note that the proposed filter (b) and the filter (c) use the same decompressed 
data iii and i/g for measurements. However, in the proposed filter, the compression error information (13) 
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and (15) are explicitly used in (18), whereas they are ignored in the filter (c). In summary, the estimation 
error reduction was possible because of two facts: (1) the proposed compression method provides the 
compression error bound (13) and (15), and (2) the proposed filter algorithm explicitly uses the error 
compression bound. 

6. Conclusion 

In this paper, we have proposed a compression method for inertial and magnetic sensor data. The 
proposed compression method guarantees that the compression error is bounded by a prescribed Cmax 
value. Smaller e^ax value usually increases the bit rate and reduces the estimation error of the filter when 
the decompressed data is used. Thus Cmax plays the role of a trade-off parameter between the bit rate 
and the estimation error. Also we have seen that by using a bound on compression error, the estimation 
error can be reduced. 
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